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Abstract 

I report on an implementation of an algorithm for the automated numerical calculation of 
spin- and colour-conelated Bom matrix elements in QCD. These spin- and colour-con^elated 
matrix elements are needed for NLO calculations in combination with the subtraction method. 
Both massless and massive quarks are considered. There are no restrictions on the number of 
external particles. As a trivial sub-case, the algorithm also applies to Born matrix elements 
without any coiTclations. These are sufficient for leading order calculations. 



1 Introduction 



QCD processes will constitute the bulk of events at the LHC. These processes provide informa- 
tion on the strong interaction and form quite often important background for searches of new 
physics. An accurate description of jet physics is therefore mandatory. Although jet observables 
can rather easily be modelled at leading order (LO) in perturbation theory [1-11], this description 
suffers several drawbacks. A leading order calculation depends strongly on the renormalization 
scale and can therefore give only an order-of-magnitude-estimate on absolute rates. Secondly, at 
leading order a jet is modelled by a single parton. This is a very crude approximation and over- 
simplifies inter- and intra-jet correlations. The situation is improved by including higher-order 
corrections in perturbation theory. 

At present, there are many next-to-leading order (NLO) calculation for 2 ^ 2 processes at 
hadron colliders, but only a few for 2^3 processes. Fully differential numerical programs 
exist for example for pp 3 jets [12-14], pp —>■ V + 2 jets [15], pp — > tiH [16, 17] and pp 
H + 2 jets [18,19]. 

It is desirable to have NLO calculations for 2 — > n processes in hadron-hadron collisions 
with n in the range of /i = 3, 4, ...,6,7. QCD processes like — > n jets form often impor- 
tant backgrounds for the searches of signals of new physics. However, the complexity of the 
calculation increases with the number of final state particles. To overcome the computational 
limitation, there have been in the past years several proposals for the automated computation 
of next-to-leading order observables [20-34]. These publications focussed mainly on the auto- 
mated computation of loop integrals. Equally important is the computation of the real emission 
contribution. It is well known that in general an NLO observable will receive contributions from 
the virtual corrections and the real emission part. Taken separately, each of the two contribu- 
tions is divergent due to the presence of infrared singularities. Only the sum of the two is finite. 
There are several general method available to handle this problem, like the phase-space slicing 
method [35-37] or the subtraction method [38^3]. In this paper I will focus on the dipole sub- 
traction method [39-43]. The dipole subtraction method requires the calculation of spin- and 
colour-correlated Born matrix elements. In this paper I describe a method for the automated 
calculation of these quantities. While the kinematical part of the matrix elements is calculated 
numerically, colour-correlation matrices are calculated symbolically at the initialisation phase of 
the program. The C++ library "GiNaC" [44] allows to mix numerical and symbolical code in a 
single program. The program uses standard techniques like spinor methods [45-49] and colour 
decomposition [50-56]. The program computes helicity amplitudes, which are decomposed into 
colour factors and partial amplitudes. The partial amplitudes are computed with the help of 
Berends-Giele type recurrence relations [1, 57]. It should be noted that recently interesting new 
methods emerged for the computation of partial amplitudes [58-61]. 

This paper is organised as follows: In the following section I present the general setup for the 
dipole subtraction method and review a few basic tools for the calculation of QCD amplitudes. 
Sect. 3 describes the algorithm for the calculation of colour-correlated Born matrix elements. 
The numerical implementation is discussed in sect. 4. Finally, sect. 5 contains the conclusions 
and an outlook. In an appendix I summarise the colour-ordered Feynman rules and the colour- 
correlation operators. Furthermore, I give some technical details on the implementation into a 
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C++ program. 



2 General setup and basic tools 
2.1 The dipole formalism 

The starting point for the calculation of an infrared safe observable O in hadron-hadron collisions 
is the following formula: 

(O) = /rfx,/(.,)/*2/fe)^p^pjj^-i- (1) 



/ 



(f) (271 + 1) {2J2 + l)nin2 

d^n{pi,P2;P3,---,Pn+2)0{pi,...,p„+2)\^„+2\ ■ (2) 



This equation gives the contribution from the n-parton final state. The two incoming particles are 
labelled pi and p2, while pi, to pn+2 denote the final state particles. f{x) gives the probability 
of finding a parton a with momentum fraction x inside the parent hadron h. A sum over all 
possible partons a is understood implicitly. 2K(s) is the flux factor, 1 / (27i + 1) and 1 / (2/2 + 1) 
correspond to an averaging over the initial helicities and ni and n2 are the number of colour 
degrees of the initial state particles, d^n is the phase space measure for n final state particles, 
including (if appropriate) the identical particle factors. The matrix element |^n+2p is calculated 
perturbatively. 

At NLO one has the following contributions: 



{Of^^ = I On+ido^ + J Onda^ + J Onda^. 



(3) 



n+1 



Here I used a rather condensed notation, da^ denotes the real emission contribution, whose ma- 
trix element is given by the square of the Born amplitudes with (n + 3) partons | *^^3 do^ gives 
the virtual contribution, whose matrix element is given by the interference term of the one-loop 

amplitude -!?^^2 ^i^-h (« + 2) partons with the corresponding Born amplitude Aj^^2- d'^^ denotes 
a coUinear subtraction term, which subtracts the initial-state coUinear singularities. Taken sepa- 
rately, the individual contributions are divergent and only their sum is finite. In order to render 
the individual contributions finite, such that the phase space integrations can be performed by 
Monte Carlo methods, one adds and subtracts a suitable chosen piece [39^3]: 



{Of^^ = J (On+ld(^ - OndCS^^ + j OndO^ + Ond<f +0nJdO^\. 

n+1 n \ 1 / 



(4) 



The matrix element corresponding to the approximation term d<s^ is given as a sum over dipoles: 

pairs i,j k^ij 
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Each dipole contribution has the following form: 

^ij,k = - X— — -^„+2 [Ph-, P[ii) ,-,Pk,-) -^^Vij,kK+2 {Ph-, P{ij) ,-,Pk,-) -(6) 

Here T, denotes the colour charge operator [39] for parton i and Vij± is a matrix in the spin 
space of the emitter parton {ij). Explicit formulae for the expressions Vij^k can be found in the 
literature [39^3] and are not repeated here. In the numerical program both the dipole terms for 
massless and massive partons are implemented. 

In general, the operators T/ lead to colour correlations, while the Vi;,/t's lead to spin correla- 
tions. The colour charge operators T/ for a quark, gluon and antiquark in the final state are 

quark: {...qi...) {Ttj) R {...qj..) , 

gluon: J?* ^(...g"...), 

antiquark : {...qi...) {-Tjl) R {...qj-) ■ (7) 

The corresponding colour charge operators for a quark, gluon and antiquark in the initial state 
are 

quark: {...qi...) {-Tf;) R {...qj..) , 

gluon: J?* (.../...) (//-^)^ (.../...), 

antiquark: il* {...q,...) [T^^ {...qj...) . (8) 

In the amplitude an incoming quark is denoted as an outgoing antiquark and vice versa. 

The subtraction term can be integrated over the unresolved one-parton phase space. Due 
to this integration, all spin-correlations average out, but colour correlations still remain. In a 
compact notation, the result of this integration is often written as 

dc5'^ + Jdo^ = l®do^ + K®do^ + F®da^. (9) 
1 

The notation (g) indicates that colour correlation still remain. The term I (g) do^ lives on the phase 
space of the n-parton configuration and has the appropriate singularity structure to cancel the 
infrared divergences coming from the one-loop amplitude. Therefore da^ + 1 (8) do^ is infrared 
finite. 

The purpose of the paper is to set up a numerical program for the automated computation of 
the terms 

j (On+ida^-Onda^^ (10) 

n+l 

and 

yo„(l(g)Jo^ + K(g)Jo^+P(g)Jo^). (11) 

n 
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This requires the computation of the matrix elements with (n + 3) partons with no spin- or colour- 
correlations (implicit in do^) as well as the computation of matrix elements with {n + 2) partons 
with spin- and colour-correlations. The subtraction terms in eq. (10) involve spin and colour cor- 
relations. The insertion operators I, K and P induce colour correlations, but no spin correlations. 
One is therefore naturally lead to the calculation of colour-ordered amplitudes in a helicity basis. 
Basic techniques for such a task are reviewed in the next subsections. 

2.2 Double line notation 

In QCD one deals with quarks and gluons. Both types of partons carry information on the colour 
degrees of freedoms and the kinematical degrees of freedom. Quarks have a colour index i, 
running from 1 to and corresponding to the fundamental representation of SU (N). The kine- 
matical information can be represented for massless quarks by Weyl spinors pA or pg, where the 
indices A or 5 run from 1 to 2. The corresponding information for gluons is in the conventional 
approach represented by a colour index a, running from 1 to A'^ — 1 and which corresponds to 
the adjoint representation of SU (N). The kinematical information is represented by a Lorentz 
index //, running from to 3. It is useful, to treat quarks and gluons on the same footing. To this 
aim, I follow the "double-line"-approach [62] and convert a gluon index to two quark indices. I 
do this for the colour degrees of freedom, as well as for the kinematical parts. 

In detail, this is done as follows: In Feynman diagrams one distinguishes edges and vertices. 
Edges are propagators as well as polarisation vectors or spinors for external particles. Vertices 
are all interaction vertices. For vector-like couplings one can write 

V,E^ = V,^^E, = V, (^<,a^^^) = {-^y.<B) ^'""^v) , (12) 

which allows us to replace a contraction over // by two contractions over A and B. One can apply 
the same trick to the colour algebra: 

ya^a ^ ya^ab^b ^ ya (^^Tf^ = (^V^T^-^V"^ (^IT^iE^^ . (13) 

Again, this equation allows us to replace a contraction over an adjoint index a by two contractions 
over indices i and j in the fundamental representation. The Feynman rules for QCD in the double 
line notation are listed in appendix A. 

2.3 Colour decomposition 

In this paper I use the normalisation 

Trr^r^ = is^^ (14) 
2 

for the colour matrices. Amplitudes in QCD may be decomposed into group-theoretical fac- 
tors (carrying the colour structures) multiplied by kinematic functions called partial amplitudes 
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[50-54]. These partial amplitudes do not contain any colour information and are gauge-invariant 
objects. 



The colour decomposition is obtained by replacing the structure constants f^^^ by 



■^abc ^ 2 



Tr ^T'^T^T''^^ — Tr ^7'^7'''T''^ 



(15) 



which follows from [7", 7^] = if^^^T^. The resulting traces and strings of colour matrices can 
be further simplified with the help of the Fierz identity : 

TtjT& = ^(sAvt-^S,A/)- (16) 
In the pure gluonic case tree level amplitudes with n external gluons may be written in the form 

VV^/ aes„/z„ 

where the sum is over all non-cyclic permutations of the external gluon legs. The quantities 
A„(oi, ...,o„), called the partial amplitudes, contain the kinematic information. They are colour- 
ordered, e.g. only diagrams with a particular cyclic ordering of the gluon s contribute. The choice 
of the basis for the colour structures is not unique, and several proposals for bases can be found 
in the literature [55,56]. Here I use the "colour-flow decomposition" [56]. As a further example 
I give the the colour decomposition for a tree amplitude with a pair of quarks: 

Jln+2{q,l,2,...,n,q) = (^-^ J^6,y„^5;„j^-„^...5,„j-A„+2(^,Oi,02, ...,0„,^). (18) 

where the sum is over all permutations of the gluon legs. In squaring these amplitudes a colour 
projector 



(19) 



has to applied to each gluon. 

While the colour structure of the examples quoted above is rather simple, the colour decom- 
position can be become rather involved for amplitudes with many pairs of quarks. A systematic 
algorithm for the colour decomposition and the diagrams contributing to a single colour structure 
is given in sect. 3. 



2.4 Spinor techniques 

For the calculation of helicity amplitudes [45^9] one chooses for the spinors corresponding 
to external massless quarks two-component Weyl spinors. Two notations for Weyl spinors are 
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commonly used in the literature. The relation between the bra-ket notation and the notation using 
dotted and undotted indices is as follows: 

\p+)=Pb, {p+\=Pa, (20) 
b-)=/, {P-\=P^- (21) 

Spinor products are denoted as follows: 

{pq) = {P- \q+) = P^qA, [qp] = {q+ \p-) = qAP^- (22) 
For the polarisation vectors of the external gluons one uses 

where k is the momentum of the gluon and q is an arbitrary light-like reference momentum. In 
the "double-line" notation this becomes 

4''(^,^) = t4v^/. e^'iKq) = -^yk'. (24) 
{qk) [kq] 

For spinors corresponding to massive quarks the formulae from ref. [63] are used. 

2.5 Recurrence relations 

Recursive techniques [1,57] build partial amplitudes from smaller building blocks, usually called 
colour-ordered off-shell currents. Off-shell currents are objects with n on-shell leg and one addi- 
tional leg off-shell. Momentum conservation is satisfied. It should be noted that off-shell currents 
are not gauge-invariant objects. Recurrence relations relate off-shell currents with n legs to off- 
shell currents with fewer legs. 
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For the pure gluon current J^^, the recurrence relation reads 



xyDCFEHGiPlJ,Pj+l,n)P'^"'^\pi,n) 

+ E E [pi,---,pf'^h---,qj)Jk-j[P%v---,Pk'^j+i^---^^k) 
j=ik=j+i ^ ^ ^ ^ 

'<Jn"k {pt+v-.Phqk+X.-.qn) VDCFEHGJiP"^''iPl,n) (25) 

This relation is pictorially shown in fig. 1. In this formula, the qt's are the reference momenta 
for the external gluons, p^^^^ (^) is the expression for the gluon propagator and ygADCFE (^i ' ^2) 
^BADCFEHG expressions for the three-gluon and four-gluon vertices, respectively. I 

further used the notation 

Phi = T,Pi- ^^^^ 

i=i 

The recursion starts with the current with one external leg, which is given by the polarisation 
vector: 

Jf{pf;qi) = 4^(pi,?i) (27) 

Similar recurrence relations can be written down for the quark- and antiquark currents, as 
well as the gluon currents in full QCD. The guiding principle is to follow the off-shell leg into 
the "blob", representing the sum of all diagrams, and to sum on the r.h.s of the recurrence relation 
over all vertices involving this off-shell leg and off-shell currents with less external legs. 



3 The method 

In this section I describe in detail the method for the automated computation of Bom matrix 
elements in QCD The matrix elements may or may not involve spin and/or colour correlations. 



3.1 Helicity amplitudes and spin correlations 

The program computes helicity amplitudes. For a given set of external momenta, each helicity 
amplitude evaluates to a complex number. If no spin correlations are present, the matrix element 
is simply given as the squared modulus of the amplitude summed over all helicity configurations. 
In the dipole formalism, spin correlations are related to the splittings g — > gg and g — > qq. In the 
original formulation of Catani and Seymour they are written as 

si;{...,p^ij),...)S^Jl4...,P(^ij),...), (28) 
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where Jl/^ denotes the amplitude with the polarisation vector of the emitter gluon {ij) amputated. 
Furthermore, the spin correlation tensor is of the form 



and the vector satisfies 

v-p^ij) = 0. 

Within the helicity formalism the spin correlation is evaluated as [64] 



(29) 



(30) 



(31) 



where Jl{...,p^^jy...) denotes the helicity amplitude, where the emitter gluon has "+", respec- 
tively "— " helicity. E is given by 



In eq. (32) q is as usual an arbitrary null reference momentum. 



(32) 



3.2 Amplitudes with more than one quark-antiquark pair 

If more than one quark-antiquark pair is present, we have to sum over all quark permutations. 
An amplitude with Uq quark-antiquark pairs can be written as 

A {qi,qi,...,q2,q2,...,qng,qng) = 

E (-l)"" f n^Sl) I ^ (^hqa{l),-,^2,qa{2),-,qn,,qo{n,)) ■ (33) 
oe5K) \j=l 7 

Here, (— 1)*^ equals —1 whenever the permutation is odd and equals +1 if the permutation is 
even. In SI each external quark-antiquark pair {qj,qc{j)) is connected by a continuous fermion 

line. The flavour factor 6{'"^, , ensures that this combination is only taken into account, if qi and 
q^si^j) have the same flavour. 

3.3 The colour structure 

The amplitude k is decomposed into colour factors and partial amplitudes: 

k = Y,CiAi (34) 

i 

Each partial amplitude A; has a fixed cyclic ordering of the external legs. For Bom graphs we 
can take this ordering such that a quark follows immediately its corresponding antiquark in the 



9 



k + 3 



I 




n 



3 



2 



1 



k 



Figure 2: The cyclic order of a partial amplitude. Without loss of generality we can assume that 
quarks follow immediately antiquarks in the clockwise order. 



clockwise orientation. This is shown in fig. 2. That is to say, that gluon are emitted from a quark 
line only to the right when following the fermion line arrow. If a gluon would be emitted to the 
left, we could draw an equivalent diagram by flipping the off-shell current attached to this gluon 
to the right of the fermion line. 

All possible cyclic orderings are generated as follows: We assume that the amplitude has Ug 
external gluons, Uq external quarks and therefore necessarily also Uq external antiquarks. We first 
note, since a quark follows immediately its corresponding antiquark, we can treat an adjacent 
((^,(5')-pair as an external "pseudo-leg", which is permutated together. The amplitude has there- 
fore Hg + Hq pseudo-legs. Then aU possible cyclic orderings are obtained by summing over all 
permutations of the pseudo-legs and factoring out the cyclic permutations, e.g. each ordering 
corresponds to an element of 



This is equivalent to fixing the first external pseudo-leg and summing over all permutation of the 
remaining (n^ + — 1) external pseudo-legs. Therefore there are 



inequivalent cyclic orderings. 

For the pure gluon amplitude {Uq = 0) each cyclic ordering corresponds to one colour factor 
c/. The situation is different if quarks are present (uq ^ 0). This is related to the fact that the 
gluon propagator in an SV (N) gauge theory can be written as a propagator corresponding to an 
U (N) gauge theory minus a part which subtracts out the additional U{1) piece. The kinematic 



S{ng + nq)/Z{ng + nq). 



(35) 



{ng + nq-l)\ 



(36) 
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UiN) 

[/(I) 




UiN) 



Uil) 



Figure 3: An example for the decomposition into colour clusters. 



parts of the U (n) and the f/(l) pieces are the same: 



jBADC 



However, they differ by their colour structure: 



UiN): 
Uil): 



I 

1? 



1 



(37) 



(38) 



Note that each propagation of a i7(l) gluon is accompanied by a factor (— 1)/A^. It can be shown 
that the Uil) gluon couples only to quark lines [56]. Therefore an amplitude with Uq quarks 
can contain up to (n^ — 1) gluons of type f/(l). Each f/(l) gluon separates a Born amplitude 
into colour-disconnected pieces. We define a colour cluster as a part of an amplitude, which 
is connected to the rest of the amplitude only by an Uil) gluon and which does not contain 
by itself any t/(l) gluon. This concept is illustrated in fig. 3, which shows a diagram with 
three quark-antiquark pairs, one U iN) gluon and one Uil) gluon. This diagram has two colour 
clusters, formed by the particles (1,2,3,4) and (5,6), and separated by the f/(l) gluon. For 
an amplitude with Uq quark-antiquark pairs one can have from 1 to Uq colour clusters. From 
the cluster decomposition the colour structure can be read off easily. The example in fig. 3 
contributes to the colour structure 



'273 



) [K.) 



(39) 



In general, given a colour cluster assignment, the corresponding colour factor c/ is constructed 
as follows: First of all, the colour factor factorizes into a product of the contributions from the 
individual colour clusters. 



Y \ {^cluster 1 ) ^cluster 

X n ^'j- 



(40) 
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Cij is the colour factor corresponding to cluster j. For a cluster consisting only of gluons, Cij is 
given by 

ghgi, •••,gn : Cij = ^hhKh-K-iAh- (41) 

An antiquark-quark pair can be treated effectively as a single gluon. For example the colour 
factor associated to a colour cluster consisting of an quark-antiquark pair and {n — 2) gluons is 
given by 

qhq2,g3, ■■■,gn ■■ Cij = ^i2j3^i3j4--\-iAh ■ (42) 
As a further example we quote the colour factor for a cluster with two quark-antiquark pairs: 

quq2,g3, ...,qk,qk+h ■■■,gn ■ Cij = ^hhh^u--\_ih\+xjk+2---K-\Ah- (43) 

The pattern should be clear. The colour factor associated to a individual colour cluster is just 
a sequence of Kronecker 5's, corresponding to the cyclic ordering of the legs belonging to this 
colour cluster. 

It remains to derive a method, how all possible colour clusterings can be generated. This 
is a combinatorial problem. For a fixed cyclic ordering we can generate all possible colour 
clusterings as follows: We first sum over the number of possible colour clusters. Let n^-iuster be 
the number of colour clusters, where Uduster ranges from 1 to Uq. For a fixed Uduster we then sum 
over all partitions of (n^ + n^) into riduster pieces rfj^^^^, such that 

^ clust er 

£ nf'^' = ng + nq. (44) 

7=1 

For a partition we take into account the order, such that for example (1,1,2), (1,2,1) and (2,1,1) 
are distinct partitions of 4. rf!-'^^^^ gives the number of external pseudo-legs belonging to cluster 
j. Obviously, an adjacent antiquark-quark pair has to belong to the same colour cluster, therefore 
it is counted as one external pseudo-leg. Finally, we have to sum over all possible starting points 
of the colour clusters with respect to the cyclic ordering. Here we observe that the members of 
a colour cluster need not be adjacent in the cyclic ordering. An example for a colour assignment 
in the cyclic ordering would be 

{q, qhg, g; g, {q,q) , g; §, gj g, (q, q) ,g,gy (45) 

clusterl cluster! clusterl clusters 

In this example, cluster 2 is embedded in cluster 1. The summation over the starting points has 
to full-fill the following requirements: 

(i) The external pseudo-leg 1 belongs to colour cluster 1. 

(ii) The colour cluster {j +1) starts after colour cluster j for all j > 2. (Colour cluster 1 may 

start at the end of the cyclic ordering.) 
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(iii) If the assignment of external pseudo-legs to colour cluster j has been interrupted by the 
starting of a new cluster k ( with ^ > 7 ), the assignment to cluster j cannot be continued 
until all members of cluster k have been assigned. 

Requirement (iii) ensures that we cannot have a sequence like cluster 1, cluster 2, cluster 1, 
cluster 2. The assignment of the external pseudo-legs to colour clusters is now done as follows: 
Let 

{mr\mr'-\....mZlJ (46) 
be an {nduster - 1) -tuple, such that 

^start < ^start ^4^^ 

and 

1 < mf' < 2 - + £ nf"'''' (48) 
k=i 

Then 



^start ^ ^start^j_^ (49) 

defines the starting point of cluster j for j = 2, nduster- The starting points nj"''' together with 
the rules (i) and (iii) define uniquely the assignment of the external pseudo-legs to the colour 
clusters. Summing over all {nduster — 1) -tuples in eq. (46) subject to the constraints (47) and 
(48) generates all possibilities with nduster colour clusters, in which colour cluster j has rfj^^^^ 
external pseudo-legs. 

Since each colour cluster couples to the rest of the amplitude through a t/(l)-gluon, it has 
to contain at least one quark- antiquark pair. Therefore configurations, where a colour cluster 
does not contain a quark-antiquark pair are vetoed, with the trivial exception of the pure gluon 
amplitude, which consists of one colour cluster and no quark-antiquark pairs. 



With the colour cluster decomposition and a method for the generation of all cluster decom- 
position at hand, 1 now turn back to the computation of the amplitude squared. From eq. (33) 
and eq. (34) it is clear that we can write any amplitude in the form 

= I^cA-, (50) 

where the c/s are the colour factors and the A,'s are the partial amplitudes which contain the 
kinematical information. In squaring the amplitude we obtain 

= Y^Ai{ciPc^^A). (51) 

hi 
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The colour projector is given as a product with one factor for each external particle: 



rig+lUq 

P = Yl Pk, (52) 

k=l 

where the individual colour projectors for a quark, antiquark and a gluon are: 

Pq = h. Pq = ^jp Ps = ^i^jj-]^lf>ji- (53) 

The only non-trivial piece is given by the colour projector for the external gluons, which is a 
consequence of the double-line notation. Note that 



(c.-Pcj) (54) 



defines a matrix, which is independent of the four-momenta of the particles. Therefore this 
matrix can be calculated at the initialisation phase of the program. As each entry is given as a 
contraction of Kronecker 5's, this can be done easily symbolically with the rules 

^ij^jk = 5*, 5,7 = A^. (55) 

The program uses the C-i-i- library "GiNaC" for this task. In appendix (C.l) I give a small 
example program. The resulting expression is a function of A^, and after substituting N = 3 the 
result can be converted to a double precision number. Note that run-time performance is not an 
issue here, since this calculation occurs only at the initialisation phase of the program. To obtain 
the amplitude squared, the matrix Mij is first calculated at the initialisation phase and stored in 
memory. Then for each momentum configuration the vector of partial amplitudes A = (A i , A2 , . . . ) 
is computed. The amplitude squared is then given by 

I^P = AMA^ (56) 

The inclusion of colour-correlations is rather straightforward. To include colour-correlation be- 
tween particles a and b, one replaces Pa and Pt in eq. (52) by the appropriate colour-correlation 
operator. For example, the colour-correlation operator • for a quark-antiquark pair reads 



n 



-n 



72 > ^ ' — 72 

A complete list of all relevant colour-correlation operators can be found in the appendix (B). The 
corresponding matrices Mij depend now on a and b, but are still independent of the four-momenta 
of the particles. Therefore they can be computed at the initialisation phase of the program. For a 
matrix element with n = ng + luq external particles, there are 

\n{n-l) (58) 

possibilities of choosing the colour-correlated partons a and b. Therefore the initialisation phase 
of the program computes and stores n{n—l)/2 different colour matrices Mij . For realistic values 
of n, say n <9, the CPU time and memory requirements for this task are rather modest. 
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3.4 The partial amplitudes 

It remains to discuss how the partial amplitudes A/, entering eq. (50) and eq. (51) are com- 
puted. This is done with the help of off-shell currents and recurrence relations. Compared to the 
introductory discussion of the pure gluonic off-shell current in section 2.5 there are additional 
complications: First of all, a rather trivial extension is given by the fact, that in full QCD we have 
to allow for the possibility of multiple quark-antiquark pairs. Secondly, and more important is 
the fact that the recurrence relations have to respect the decomposition of a partial amplitude into 
colour clusters. The algorithm is summarised as follows: 

(i) We consider coupled recurrence relations for the off- shell currents corresponding to an 

i7(A'^)-gluon, an i7(l)-gluon, quarks and antiquarks. Note that the i7(A^)-gluon and the 
U ( l)-gluon are treated separately, as the latter couples only to quarks. It is also convenient 
to distinguish the off-shell currents for the quarks and antiquarks, depending on whether 
the quarks are massive or massless. In the latter case specialised (and faster) routines can 
be used, since only helicity conserving interactions enter. 

(ii) All recurrence relation express an off-shell current of type A as a sum over off-shell currents 

with fewer legs, which are combined through the basic three- or four-valent vertices of the 
theory. The recurrence relations takes into account all possible interaction vertices, which 
contain A. Note that the off-shell currents, which enter the r.h.s. of the recurrence relation 
need not be of type A. For the example, the recurrence relation for an U (A'^)-gluon involves 
the quark-antiquark-gluon vertex and therefore the off-shell currents for a quark and an 
antiquark. In general, the recurrence relations yield a coupled system of equations. 

(iii) The off-shell parton for the quark-current, the antiquark-current and the ?7(A'^)-gluon- 
current belongs to a specific colour cluster a. The recurrence relation splits the off-shell 
current with n external legs into off-shell currents with less external legs. This splitting has 
to respect the following selection rules: 

- For the U (A^) -current, the off-shell current attached through the three- and four-gluon 

vertex have to contain at least one parton belonging to colour cluster a. In the off-shell 
quark- and antiquark-current, which are attached through the gluon-quark-antiquark 
vertex to the U (A^)-gluon current, the off-shell quark- and antiquark-lines have to 
belong to colour cluster a. 

- For the off-shell quark current, the sub-current attached through an U (A^)-gluon must 

contain at least one parton belonging to colour cluster a. On the other hand, the sub- 
current attached through an i7(l)-gluon may not contain any parton of colour cluster 
a. Similar considerations apply to the antiquark current. 

- Finally, the off-shell i7(l)-current is rather simple and the recurrence relation involves 

an quark- and an antiquark-current, whose off-shell legs necessarily belong to the 
same colour cluster. 

(iv) As a further selection rule we have to veto configuration, where the off-shell current is 
divided into sub-currents between leg j and (7 + 1), in the case where these two legs 
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belong to the same colour cluster b, which is different from the colour cluster a of the 
off-shell leg. That is to say, that the recurrence relation where the off-shell leg belong to 
cluster a, cannot split legs which belong to a different colour cluster b. 

3.5 The pure gluon amplitude 

In principle, the pure gluon amplitude can be treated with the methods discussed above. However 
the pure gluon amplitude is a rather special case, which leads to many additional simplifications. 

Since it is known that pure gluonic processes will contribute significantly to the cross section at 
the LHC, it is desirable to treat these processes separately with optimised routines, taking into 
account the additional simplifications. The simplifications are: 

- There is only one colour cluster and the colour decomposition is simply given by the {ug — 1) ! 

inequivalent cyclic orderings, as in eq.(17). 

- ?7(l)-gluons can be ignored and the recurrence relation for the partial amplitudes is given by 

eq. (25). 

- In calculating the colour matrix Mjj, the colour projectors Pg in eq. (53) may be replaced by 

Pg ^ dj^djj. (59) 

3.6 QCD amplitudes with one electro-weak boson 

The methods discussed above require only minor modifications to include amplitudes with QCD 
partons and one electro-weak boson. As these are relevant to electron-positron annihilation, 
electron-proton collisions or Z-production at the LHC / Tevatron, these amplitudes have been 
implemented as well. The amplitudes are computed by considering a recurrence relation, which 
couples the electro- weak current to an off-shell quark current and an off-shell antiquark current. 

4 Numerical implementation 

The algorithms discussed above have been implemented into a computer program. This numer- 
ical program can compute Bom matrix elements in QCD with spin- and colour-correlations. To 
test the program I have first considered the case, where spin- and colour correlations are absent. 
In this case one can compare the results with the ones from the program Madgraph. 1 quote here 
the results of this comparison for processes with up to seven external particles. The labelling of 
the momenta is 

PlP2^ P3,P4,-,Pn- (60) 

pi and p2 are the incoming momenta, p3 to pn are the outgoing momenta. For 2 — > 2 processes 
I took the following set of momenta (in units of GeV): 

pi = (45.0,0.0,0.0,-45.0), 
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P2 = (45.0,0.0,0.0,45.0), 

p^ = (45.0,-20.8997,-29.6778,26.5976), 

P4 = (45.0,20.8997,29.6778,-26.5976). (61) 

The same initial state momenta pi and p2 are used for all other processes. The final state 
momenta for the 2 — > 3 processes were chosen as 

P3 = (41.8145,-9.20663,-26.7503,30.7914), 

P4 = (17.3829,12.8067,10.7712,-4.70487), 

P5 = (30.8026,-3.6001,15.9791,-26.0865). (62) 

For 2 — > 4 processes I used 

P3 = (29.9152,-18.1846,-8.69254,22.1061), 
P4 = (9.82719,4.07529,8.79524,-1.61538), 
P5 = (22.171,-9.26417,14.187,-14.2988), 

P6 = (28.0866,23.3735,-14.2898,-6.19197). (63) 

Finally, for 2 ^ 5 processes I used 

P3 = (20.165,-13.0392,0.0298292,15.3819), 

P4 = (9.60811,2.4114,9.15728,-1.6264), 

P5 = (20.5589,-7.64505,15.4771,-11.166), 

P6 = (18.087,17.056,-3.25968,-5.06046), 

pj = (21.581,1.21688,-21.4045,2.47093). (64) 

The strong coupling constant was taken to be = 0. 1 18. For this comparison, all quark masses 
have been set to zero. The flavour labels serve only to distinguish identical quarks from non- 
identical quarks. Table 1 shows the comparison of our program with Madgraph for the computa- 
tion of the matrix elements corresponding to the indicated processes. The results do not contain 
any averaging over the colour degrees of freedom for the initial-state particles, nor do they con- 
tain symmetry factors for the final-state particles. As can be seen from the table, the agreement 
is satisfactory. 

To check spin- and colour-correlations 1 have compared the program with existing NLO codes 
for e+e^ 4 jets [64] and pp tig [65]. 

Table 2 gives an indication for the CPU time needed to evaluate matrix elements of increasing 
complexity. It gives the CPU time needed for the computation of the matrix elements, summed 
over all colours and spins, corresponding to the following cases: The amplitude J? (gi , g„) with 
n gluons, the amplitude ^{q,q,g^, ...,gn) with an q, q-^dix and {n — 2) gluons and the amplitude 
Sl(q,q,q',q',g5,..., gn) with two distinct q, ^-pairs and {n — 4) gluons. 

5 Conclusions and outlook 

In this paper I discussed an algorithm for the automated computation of spin- and colour-correlated 
Bom matrix elements in QCD. These matrix elements are needed for NLO calculations in com- 
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Process 



gd_-^d_g 
ud du 

dd ^ dd 



gd_^dgg 
ud — > dug 
dd — > ddg 



gd_^d_ggg 
ud — > dugg 
dd — > ddgg 
ud ^ duss 
ud — > duiiu 
dd — > dddd 



ud — > dusgs 
ud duugu 



this work 



56203.4 
8436.64 
1374.01 
1287.74 



21269.2 
3222.01 
56.459 
53.2424 



1354.24 
138.691 
0.975563 
0.902231 
0.0116469 
0.0524928 
0.0583822 



0.000453678 
0.00202449 



Madgraph 



56203.2 
8436.62 
1374.01 
1287.74 



21269.3 
3222.02 
56.4591 
53.2425 



1354.22 
138.689 
0.975546 
0.902215 
0.0116467 
0.0524927 
0.0583821 



0.000453671 
0.00202446 



Table 1: Comparison of our program with Madgraph for various matrix elements with up to 
seven external particles. 



n 


4 5 6 7 8 


time for 
time for 
time for 


^{g\,-,gn)? 
Jl{q,q,q',q',g5,...,gn)\^ 


0.0006 0.009 0.18 4 127 
0.0004 0.003 0.05 0.6 14 
0.0002 0.002 0.02 0.4 8 



Table 2: CPU time in seconds for the computation of some matrix elements summed over all 
helicities and colours on a standard PC (Pentium IV with 2 GHz). The examples consists of the 
n gluon amplitudes, the amplitudes with an q, q-pmr and (n — 2) gluons and the amplitudes with 
two distinct ^,^-pairs and {n — 4) gluons. 
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bination with the subtraction method. I implemented the algorithm into a computer program. 
The program handles QCD amplitudes with massless and/or massive quarks. In addition, I have 
implemented the extension to QCD amplitudes with one additional electro-weak boson. 

The methods presented here are part of a larger project for the automated computation of 
observables at next-to-leading order for LHC physics. The remaining missing piece is the auto- 
mated computation of the interference term of the one-loop amplitude with the Bom amplitude. 
In a previous publication, we already reported on the automated computation of the one-loop 
integrals entering the one-loop amplitude [31]. Work on the automated computation of the inter- 
ference term is in progress. 
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A Feynman rules 

In this appendix I summarise the colour-ordered Feynman rules. I extract from each formula the 
coupling constant and split the remainder into a colour part and a kinematical part. 



A.l Propagators, polarisation vectors and polarisation sums 
Gluon propagator 

In Feynman gauge, the gluon propagator is given by —ig^^^^^/k^. Contraction of the kinematical 
part -ig^"" jk^ with (l/2)a^^^a^'^^ yields: 

pBAbC^^^ = BA ^QQQQj DC = ^ (-E^^E'^cj (^5) 

The colour factor 5"^ is contracted within the double-line notation with ^/iT^j^/lT^f 

v/27;.^5«*v^r* = huhj-ljhjhi (66) 

The colour structure is split into two pieces. The first piece ^u^kj corresponds to the propagation 
of a U{N) gluon, whereas the second piece —hijhki/N subtracts out the additional U{1) gluon. 
Schematically we have 



i . / 

]^--^[ = 4^.A/- (67) 



Note that each propagation of a i7(l) gluon is accompanied by a factor {—\)/N. 
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Quark propagator 

The kinematical piece of the quark propagator reads: 

i 



The colour factor is simply 



(68) 



' ^ ^ = (69) 

Gluon polarisation vectors and colour projector 

The gluon polarisation vectors are given by 

4^(^,^) = A^^""' ^^'ik,q) = j^yk''. (70) 
{qk) [kq\ 

k is the momentum of the gluon and q is an arbitrary light-like reference momentum. The depen- 
dence on q drops out in gauge-invariant quantities. 

Colour factor: In the conventional approach we sum for the squared matrix element for each 
gluon over all eight colour degrees of freedom. In the double-line notation a factor v^T-" is 
moved at each end into the colour projector. Therefore, the colour projector reads 

VlTtjVlTl^i = hnhj-]^^ijhi- (71) 

A.2 Vertices 
Quark-gluon vertex 

The kinematical part of the quark-gluon vertex is given by 

CD CD 

^ I =-iV2zcAtDB^ A , I , g =-/y26/6/. (72) 

A B 

The colour factor is given by 



Ik 

1 

V2 



—^ahj. (73) 



Here I neglected terms proportional to hki, which vanish when contracted into the gluon propa- 
gator. 
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Three gluon vertex 

The kinematical part of the three-gluon vertex is given by 




AB 



The colour factor reads: 



(74) 



73 ' 3 72 '2 



(75) 



Four gluon vertex 

The kinematical part of the four-gluon vertex is given by 



GH AB 




= 2/ \lZAE^BF^CG^bH ~ ^AC^Bb^EG^FH ~ ^AG^BH^CE^bf-] ■ (^6) 



EF CD 
The colour factor reads: 



UJ4. 



(77) 



73^3 72^2 



B Colour correlations 

In this appendix I list all colour-correlation operators Tq ■ Tt, between two partons in the double 
line notation. 

j?*(...a,...,Z>,...)(T«-Ti)^(...a,...,Z>,...). (78) 

Their action between amplitudes is defined in eq. (7) and eq. (8). As we write all amplitudes in 
the colour-flow decomposition, we would like to know the action of these operators in this basis. 
In the following I denote the colour indices of the amplitude M * with barred indices, the colour 
indices of the amplitude with un-barred indices. 
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Quark-quark ( ) 



12- 



o 

i 

o 
o 



2 ( ^'i'2^i2'i j^^h'i^hh 



-12 



Quark-antiquark ( ) 



J2 



o 



o 
o 



2 (^'iJ2^J2h j^^hh^jih 



-J2 



Antiquark-antiquark (Tq-Tq) 



Ji 



J2- 



-Jl 



-J2 



2 i^jiji^jiji j^^jiji^jiJi 



Quark-gluon ( ) 




'2,72 



Antiquark-gluon ( To T„ ) 



71 



-7i 



12 , 72 12 , 72 



Gluon-gluon ( T„ T„ ) 



'1,71 



'2,72 




'1,71 



'2,72 



2 i^niAAl^jin -^'AAiAh 
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C Details on the implementation 



In this appendix I provide some details on the implementation of the algorithm into a C++ pro- 
gram. I will discuss in a small example how the colour algebra is performed. I will also give 
some hints on the implementation of loops over multi-indices like permutation, partitions, etc.. 

C.l Colour algebra 

Below I show a small program, which defines the colour structures 

ci = ^hjihh . A = ^jih . (85) 

and contracts them: 

cn = cicj. (86) 
The result is obviously c\i = N^, which equals 9 for A/^ = 3. 

#include <iostream> 

♦ include ' 'ginac/ginac.h' ' 

int main{) 
{ 

using namespace GiNaC; 

// number of colours 
int No = 3; 

// define colour indices 

ex il = idx( symbol ( "il" ) , No ) ; 

ex 12 = idx{ symbol { "12" ) , Nc ); 

ex jl = idx ( symbol ("jl") , Nc ); 
ex j2 = idx( symbol ("j2") , Nc ); 

// define colour structures 

ex cl = delta_tensor (il, j2) *delta_tensor (12, jl) ; 

ex cl_conj = delta_tensor ( j2, il) *delta_tensor ( jl, 12) ; 

// square it and contract indices 

ex cll = cl_conj * cl; 

cll = cll . simplify_indexed ( ) ; 

// convert the result to a ^Mouble'' variable 
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double c_double = real (ex_to<numeric> { cll ) ) .to_double () ; 



std::cout << ''result = ' ' << c_double << std::endl; 



return 0; 

} 



C.2 Summing over multi-indices 

The algorithms involves the summation over multi-indices. A rather simple example for a multi- 
index would be a ^-tuple (/q, i'l, ijt-i) where each entry can take values from to — 1. Other 
examples are the sum over permutations of k elements as in eq. (35) or the multi-index in eq. 
(46). To make the code readable it is desirable to write the loop as 

{ 

int N = 7; 
int k = 3; 

multi_index i_multi (N, k) ; 

for( i_multi . init ( ) ; ! i_multi . overflow () ; i_multi++) 
{ 

// can use i_multi[0], i_multi[l], etc. here 

} 

} 

and to hide the details on how the multi-index is increased into a separate class. A possible 
header file for the class multi_index could look as follows: 

class multi_index { 



public : 

multi_index (size_t N, size_t k) ; 



// functions 

multi_index & init (void) ; // initialization 

bool overflow (void) const; // returns overflow flag 

multi_index & operator++ (int); // postfix increment 

size_t operator [] (size_t i) const; // subscripting 



// member variables : 
protected : 
size_t N; 

std: : vector<size_t> v; 
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bool f lag_overf low; 

}; 



This class contains a method in it to initialise the multi-index to the first value, an operator ++ 
which increases the multi-index to the next value and method overflow, which returns true if all 
values have been run through. 
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